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The phase diagram of the Heisenberg ferromagnetic model in the presence of a 
magnetic random field (we have used bimodal distribution) of spin 5 = 1/2 (quantum 
case) and S = co (classical case) on a simple cubic lattice is studied within the 
framework of the effective- field theory in finite cluster (we have chosen N = 2 spins). 
Integrating out the part of order parameter (equation of state), we obtained an 
effective Landau expansion for the free energy written in terms of the order parameter 
^'(m). Using Maxwell construction we have obtained the phase diagram in the T — H 
plane for all interval of field. The first-order transition temperature is calculated by 
the discontinuity of the magnetization at T* {H) , on the other hand in the continuous 
transition the magnetization is null at T = Tc{H). At null temperature (T = 0) we 
have found the coexistence field He = 3.23 J that is independent of spin value. The 
transition temperature Tc{H) for the classical case {S = oo), in the T — H plane, is 
larger than the quantum case {S = 1/2). 
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I. INTRODUCTION 



Phase transitions are one of most interesting phenomena that occurs in nature. Many 
systems have phase transitions in critical regions and it is widely known that the classic Ising 
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model (and others) displays a second order temperature driven phase transition. In particu- 
lar, phase transitions and the critical behaviors of the random field Ising model (RFIM) were 
studied extensively in the last years, see jl-3| and references therein. The RFIM leads to 
a number of challenging problems in the physics of disordered systems {2-^. There are two 
basic types of disorder in spin models: i) disordered bonds (spin-glass models) and ii) site 
disorder (randomness of the applied magnetic field in the RFIM). Mean- field theory has been 
one of several techniques used to study the RFIM. Although the mean-field version of the 
RFIM is much easier, there are some open questions about the behavior of the RFIM with 
more realistic, short-rgange interactions, which still motivate experimental and theoretical 
invest igat ions [2I, IsJ. The lower critical dimension and the existence of an ordered phase in 
the three-dimensional case, have been rigorously established by mathematical proofs js, 6|. 
However, the existence of a tricritical point (TCP) for a double-5 distribution of random 
fields, in accordance with mean-field results, is still under question( see fzj, 0] and references 
therein) . 

The RFIM is revelant for the description of several physical situations, as example: 
i) for the structural phase transitions in random alloys, ii) for the phase transitions in 
commensurate charge-density-wave systems with impurity pinning and iii) in binary fiuid 
mixtures in random porous media. Random fields have been used to mimic frustration 
introduced by the disorder in interacting many body systems and for explaining several 
aspects of electronic transport in disordered insulators [91] and in systems near the metal- 
insulator transition jl^, [ll| . On the other hand, the physics of hysteresis, of the avalanche 
behavior, and of the origin of self-organized criticality 12|], has been explained by resorting 
to the analysis of the non-equilibrium behavior of suitable RFIM. There is a new class of 
problems related to the self-generated glassy behavior, which has been explained instead in 
terms of a spin model in infinitesimal random fields 13|, and more recently, the RFIM has 
been employed to describe critical behavior of amorphous magnetic systems, such as thin 

nn 

films and critical surface behavior of the amorphous semi-infinite systems p^, 115| . 

In the last years, a new effective field theory (EFT) has been used to study second-order 
phase transition of both classical and quantum spin models, and tricritical point in the 



phase diagram, which leads to useful qualitative insights for the critical behavior. The EFT 
method uses the Callen-Suzuki identities 3] as a starting point and utilizes the differential 
operator technique, developed by Honmura and Kaneyoshi 17l|. It provides a hierarchy of 
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approximations to obtain thermodynamic properties of magnetic models. One can continue 
these series of aproximations considering increasing clusters which leads to better results. 
The exact solution would be obtained by considering an infinite cluster. However, by using 
relatively small clusters that contain the topology of the lattice, one can obtain a reasonable 
description of thermodynamic properties as it will be shown below. 



Several spin models, such as Blume-Capel [l8|, random field Isi ng |l9l.l20l| . Heisenberg 



23|, Ising metamagnet 



24 
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25| and Ising with four-spin couplings 26|, |27| models, have been 



treated by using EFT. In these works, the first-order line could not be obtained due to 
the absence of an expression for the free energy. Therefore, only second-order lines and 
tricritical points were analyzed. In particular, Fittipaldi and Kaneyoshi j2^ have used 
the EFT approach to study the phase diagram of the Blume-Capel model with spin-1 on 
a two-dimensional lattice. The position of the first-order transition was obtained from the 
isotherms in the m — H plane (where m and H are the magnetization and the magnetic field, 
respectively) applying the Maxwell equal area construction. The first-order lines obtained 
in Ref. [28] are not correct, since in the limit a = —0.50 (where a = J'/ J, and J'{J) is the 
biquadratic (bilinear) coupling) at T = the exact value is D/J = —0.75, and the value 
presented in Ref. [28] was D/J = —0.50 (see Fig. 1). 

Recently, de Albuquerque, et al. [2l|, |22| have studied the phase diagram of the random 
field classical Heisenberg model (RFHM) on a simple cubic lattice. Oubelkacem, et a/. j23| 
extended the calculation to treat the quantum spin-1/2 random field Heisenberg model and 



obtained only second-order lines and tricritical points |21l-l23|. The purpose of this work is 
to discuss the complete phase diagram (entire range of the field) in the T — H plane of the 
random field Heisenberg model on a simple cubic lattice by using EFT in two-spins cluster 
(EFT-2). 

In the present work, our goal is to propose a functional for free energy, in order to obtain 
the first-order line in the phase diagram in the T — H plane for the random field classical 
and quantum spin-1/2 Heisenberg model on a simple cubic lattice. The outline of this paper 
is as follows: the model and formalism are developed in section II, and the results and 
conclusions are discussed in section HI. 
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II. MODEL AND FORMALISM 



In order to obtain the free energy, we developed the calculations to treat the phase 
diagram of the RFHM with classical {S = oo) and quantum {S = 1/2) spins. The RFHM 
is described by the following Hamiltonian: 



(1) 



where the first sum is carried out only over pairs of nearest-neighboring sites with the 
interaction J. Also 5*^^ is the z-component of the spin operator (vector) Sj = (Sf , Sf , S-) 



at site i. For the classical case 



we consider the normalization condition ^ (5*^^ 

fj,=x,y,z 

and for the quantum case Sj is now considered as Pauli spin operator-1/2. Hi is the random 



magnetic field that obeys the following bimodal distribution: 

V{Hi) = ^ [6{H, -H) + 6{H, + H)] 



(2) 



in which H = {Hf)^ is the root mean square deviation of the magnetic field correspondent 
to the configurational average of the probability distribution V{Hi). 

The thermal average of a general function involving spin operator components in a finite 
cluster 0({N}) can be obtained by the generalized relation of Callen and Suzuki (l6|, i.e., 

'Tr{^J0({N})e-^^4 



(0({N})) 



Tr{A.}{e-/3^-} 



(3) 



where the partial trace Trj^v} is taken over the set of N spin variables specified by a finite- 
system Hamiltonian "Hn. (■ ■ ■) indicates the canonical thermal average taken over the 
ensemble defined by the complete Hamiltonian (1). 

The Callen- Suzuki identity for a finite cluster with two-spins was derived for the first 



time by Bobak and Jascur [30] to study the criticality of the Ising model. It has a 
generalized for the description of the quantum spin- 1/2 Heisenberg ferromagne 



antiferromagnet |32|. De Sousa and Albuquerque [33[ (see also Refs. [23 



24 



so been 



3l| and 
have 



3) 



applied EFT-2 on the classical n-vector model. Latter, the EFT-2 approach was used to 
study the magnetic properties of the quantum spin-1 Heisenberg ferromagnet {34 1. More 
recently, this new EFT has been successfully used to treat second-order phase transitions of 
classical and quantum models [35|-|39|, and also to treat first-order transitions |40l-l46|. 
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Using a two-spin Hamiltonian for the finite system ^2 in the Eq. (3) (see more details 



in Refs. 



23 



24 



26|). the magnetization per spin m = {S^ + 6*2)) is found. Applying the 
differential operator technique and EFT, an approximate expression for m is obtained for all 
values of z. In particular, for the simple cubic lattice (z = 6) case, the average magnetization 
m is given by the following expression: 

4 

A(m, T,H) = J2 Mr+i{T, H)m^''+\ (4) 



m 

r=0 



where 



A(m, T, H) = [{a^ + m(3^) ■ {ay + m(3y)f G{x, y) I^^^^q , (5) 
G(x,y) = ^[G^''(x,y) + G!:«(x,y)], (6) 
nq I \ sinh(x + y±2h) 

G^{x,y) = -1^^^=^^^= (quantum case), [7) 

cosh(x + y ± 2h) + e'^^ cosh y (x — y)^ + AK'^ 

nc ( . sinh(x + y±2h) i ^ ■ ^ \ (o\ 

G,[x,y) = — — — — — — r (classical case), (8) 

cosh{x + y ±2h) + (f){K) cosh{x - y) ^ ^' 

, , 1 - C{3K) . . 

1 ( d'Mm,T,H) \ 

and 

C{x) = coth(a;) — l/x (Langevin function) (11) 

where a^ = cosh{KDf^), = smh{KD^) (/i = x,y), = ^ is the differential operator, 
K = J/ksT and h = H/kBT. The coefficients Ar{T,H), Eq. (10), are determined by 
applying the identity e"'^''^^^^G{x, y) = G{x + a, y + b), and other corresponding expressions 
that are rather lengthy to be reproduced here. 



III. RESULTS AND CONCLUSIONS 



The EFT-2 was dev eloped 
spin 



32|(see also Refs. 



23 



'or the quantum spin- 1/2 Heisenberg 3l| and classical 



26| ) ferromagnet. Therefore, the expression from Eq. (4) 



has been obtained. This new method (EFT-2) was also used to study the criticality of the 
quantum spin-1 anisotropic Heisenberg ferromagnet 3J] . It has been observed, from these 
works, that the critical temperature ksTc/ J increases with increasing spin {S) value, i.e.. 
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ksTc/J — 1.222, 3.434, and 5.030 for S = 1/2, 1, and oo, respectively. These critical be- 
havior for the dependence of with the value of the spin 5*, our results confirm the known 
results of series expansion where the values found are ksTc/ J — 0.830, 2.72, and 4.329, 
for S = 1/2, 1, and oo, respectively. For a continuous phase transition, m{T,H) decreases 
as the temperature increases and at T = Tc{H) the order parameter is null (continuously). 
Then from Eq. (4) one can locate the second-order line through the condition 

A,{T,,H) = 1, (12) 

with A^iTcH) > 0, and, additionally, the tricritical point can be located when 

A,{T,,H) = 0, (13) 

with A^{Tc,H) < 0. Depending on the range of the ratio S = H/J, we have second-order 
(0 < 5 < 6t) and first-order {6 > 6t) transitions, where (5t,Tf) is the tricritical point. One 
can note that it is not possible to calculate first-order transition lines in the basis of only 
the equation of state (4) because in this case m 7^ at the transition point. To solve this 
problem one needs to compute the free energy for the ferromagnetic (F) and paramagnetic 
(P) phases. First-order transitions then correspond to locus on the phase diagram where free 
energies are equal. Assuming that the equation of state (4) is obtained by the minimization 
of a given free energy functional like ^(m,T,i7) (i.e., |^ = 0), we can express such relation 



dm 

as 

\ m r ^ /!_ m 1 ^ 



v^(m,T,ff) = Ao(T,//) + ^i%^ 



r + 1 

r=0 



m\ (14) 



where Ao(T, H) and Ai(T, H) are arbitrary functions which turn out to be irrelevant when 
searching for the phase transition. The Eq. (14) just represents qualitatively a Landau-like 
expansion, that can not be used to obtain the thermodynamic properties, only to study the 
phase diagram of spin system. This purpose for the free en ergy functional has been recently 



applied with success to study spin systems with frustration |40l-l46l| . In the present paper, we 
use it in the random field Heisenberg model to certify the potentiality of the methodology. 
It is known that this Landau expansion for m is given by a finite serie and it is possible to 
show that Xi{T,H) > 0. Thus, we assume that this parameter Xi{T,H) is also positive in 
Eq. (14). Near the criticality (i.e., T ~ Tc, m ~ 0) we have, from the equation of state (4), 
the behavior of the magnetization given by m ^ wi^^^^i (classical critical exponent. 
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(3 = 1/2) and, consequently from Eq. (14) |^ ~ — 2 [1 — Ai(T,H)] > that corresponds 
to a minimum point (stability limit). We note that A^lT,!!) > and Ai{T,H) > 1 for all 
H < Ht (tricritical field) and T < T^. From Eq. (14), we obtain the separation point of the 
two phases F (m ^ 0) and P (m = 0), i.e., ^F{m,T,H) = ^p{0,T,H) 

r=0 + 1 



In Refs. [23|, |26|, the Eqs. (12) and (13) have been used to obtain the critical frontier which 
separates the F phase from the P phase and the tricritical point (TCP) for the classical and 
quantum cases. Simultaneously solving Eqs. (4) and (15) we obtain the second-order line 
when m = and first-order line when m ^ 0. The corresponding phase diagram in the T — H 
plane is depicted in Figure 1 for the classical and quantum spins. As a first observation, 
we note that the nature of variations of versus H reveal a common basic behavior - the 
transition temperature decreases when H/J increases, reaching the zero temperature limit 
at same value of Hd J (i.e., Hd J = 3.23). We have also observed that Tc{H) for the classical 
case is larger than the quantum case, what is accepted physically. 

In conclusion, we observe that EFT formalism allows us to study the random field Heisen- 
berg (classical and quantum) model with correlation and phase diagram in the T — H plane. 
The results by using the functional for free energy are satisfactory to calculate the first-order 
line with qualitative and, to a certain extent, quantitative confidence. We can also extend 



the presented methodology to study the magnetic properties 
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FIG. 1: Phase diagram in the T — H plane of the random field Heisenberg model on a simple cubic 
lattice for quantum (a) and classical (b) spin cases. The solid and dashed lines correspond to the 
second- and first- order phase transition respectively. The tricritical point is marked by a back 
point. 



